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Optically measuring in the photon counting regime is a recurrent challenge in mod¬ 
ern physics and a guarantee to develop weakly invasive probes. Here we investigate 
this idea on a hybrid nano-optomechanical system composed of a nanowire hybridized 
to a single Nitrogen-Vacancy (NV) defect. The vibrations of the nanoresonator grant 
a spatial degree of freedom to the quantum emitter and the photon emission event can 
now vary in space and time. We investigate how the nanomotion is encoded on the de¬ 
tected photon statistics and explore their spatio-temporal correlation properties. This 
allows a quantitative measurement of the vibrations of the nanomechanical oscillator 
at unprecedentedly low light intensities in the photon counting regime when less than 
one photon is detected per oscillation period, where standard detectors are dark-noise- 
limited. These results have implications for probing weakly interacting nanoresonators, 
for low temperature experiments and for investigating single moving markers. 


Recent developments in hybrid mechanical quantum 
systems share the common objective of creating and ex¬ 
ploring non-classical states of motion of macroscopic ob¬ 
jects These systems combine two elementary bricks 

of quantum mechanics: a mechanical oscillator and a two 
level system, in the form of superconducting qubits 12H3, 
single spins EHni, quantum dots [IMS], BEG [201 ED, 
molecules [22] or ions. Advances in this field are oriented 
towards both increasing the hybrid coupling strength and 
reducing the nanoresonator dimensions in order to max¬ 
imize the qubit sensitivity to the nanoresonator dynam¬ 
ics. Since this reduction in size concomitantly weakens 
its interaction with standard opto- or electro-mechanical 
probe fields, it is important to investigate alternative 
nanomotion readout strategies, directly based on qubit 
state measurement. Most of the interfaced qubits can be 
probed through optical or microwave fields. We analyze 
here how the emitted photon statistics convey informa¬ 
tion on the oscillator nanomotion. 

In this article we investigate these ideas on the example 
of a NV center acting as a single photon source attached 
to the vibrating extremity of a nanomechanical oscillator 
[iin]. Its nanomotion becomes encoded onto the quan¬ 
tum emitter state once immersed in a strongly confined 
pump light field, as obtained in the focus of a high numer¬ 
ical aperture objective. We investigate how the measured 
photon statistics are impacted when the position of the 
single emitter is moving in space due to the nanomotion. 
The Brownian motion of the nanoresonator is responsi¬ 
ble for a novel photon bunching signature, consequence of 
the finite size of the measurement volume. Measuring the 
temporal cross-correlation function of the photon fluxes 
collected from different locations in space permits to in¬ 
vestigate the emitter trajectory in space and time and to 
establish connections with the autocorrelation function of 
the nanoresonator position fluctuations. Finally, we illus¬ 
trate the possibility of probing the thermal noise of the 
nanoresonator at ultralow light intensities (:^100 aW), in 
the photon counting regime, when photon counting rates 
become comparable to the oscillation frequency. 


Various mechanisms responsible for position dependent 
absorption or fluorescence rates can modify the qubit 
emission statistics, as investigated in single trapped ion 
experiments [24ll26] or proposed for hybrid nanomechan- 



FIG. 1. A single quantum emitter hybridized to 
a nanomechanical oscillator. SEM (a) and CCD fluo¬ 
rescence image (b) of the hybrid system revealing the pres¬ 
ence of a bright quantum emitter at the nanowire extremity, 
(c, d): Sketch of the experimental setup (D: dichroic, BS: 
beam splitter, APD: avalanche photodiode, QPD: quadrant 
photodiode), (e): normalized fluorescence point spread func¬ 
tion (PSF) n(r) of the apparatus measured by scanning the 
NV defect in the horizontal xz plane. A dual single photon 
counter arrangement allows to verify the emitter single pho¬ 
ton source character through intensity autocorrelation mea¬ 
surements (f). (g): Calibrated displacement noise spectrum 
of the nanowire revealing the Brownian motion of the first 
eigenmodes at 300 K and atmospheric pressure [T] (inset: dif¬ 
ferential DC transmission scan along x). 
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ical systems EaiiH]. However this work focuses on fun¬ 
damental variations of the detected photon statistics re¬ 
lated to the measurement apparatus. It can thus be 
viewed as a dynamical extension of the original exper¬ 
iments of Hanbury-Brown and Twiss [29] who derived 
the spatial coherence of a light source through static cor¬ 
relation measurements from different locations in space. 
This work is of practical interest for nano-optomechanics 
at ultra-low photon fluxes, for nanosystems weakly cou¬ 
pled to light fields and those which cannot sustain large 
optical powers, such as carbon nanotubes or cryogenic en¬ 
vironments. Furthermore, we suggest that this approach 
can directly be transposed in experiments on optically 
trapped single quantum emitters [301EB or for investi¬ 
gating the diffusion properties of biological markers. 

A single photon source with a mechanical degree of 
freedom- The quantum emitter, a single NV defect hosted 
in a 50 nm diamond nanocrystal is attached to the free 
extremity of a 46 jam long SiC nanowire with a diam¬ 
eter of 200 nm, mounted at the apex of a metallic tip 
(see Fig. 1)|6]. The hybrid system is investigated with a 
confocal microscope apparatus based on high numerical 
aperture objectives (0.75 NA) (see SI). The 532 nm pump 
laser serves for both measuring the position of the hy¬ 
bridized nanoresonator using the transmitted or reflected 
beams [1] and pumping the NV defect. Its fluorescence in 
the 630 -750 nm band is detected on avalanche photode¬ 
tectors operated in the photon counting regime, featuring 
dark noise lower than 50 counts per second 10 aW). 
Fluorescence images of the NV-functionalized nanores¬ 
onator reveal the presence of a NV defect at its extrem¬ 
ity, see Fig. lb, whose single photon source character is 
verified through autocorrelation measurements (Fig. If). 
Piezo-scanning the suspended NV defect in the waist 
area allows to determine the fluorescence point-spread- 
function n(r) of the apparatus (Fig. le), featuring a mini¬ 
mum waist of Wo 380 nm. The pump laser spot size can 
also be broadened by defocussing to produce an homo¬ 
geneous illumination over the NV oscillating trajectory, 
which can be efficiently driven through electrostatic ac¬ 
tuation, see Fig. 2. In order to collect the fluorescence 
from different positions across the emitter trajectory, the 
detectors are mounted on piezo-positioners (see SI). This 
permits to explore the photon cross-correlations from dif¬ 
ferent points in space [29]. 

The Brownian motion of the nanowire is detected on a 
quadrant photodetector (see method in m), its displace¬ 
ment noise spectrum Sx[^] = f e^^^(Sx(t)Sx(t -j- T)}dT 
is shown in Fig. Ig. The first eigenmode sits 30 dB on 
top of a shot noise limited background for 1 mW injected 
power. It oscillates at a frequency = 190 kHz 

with a quality factor of Q 2, drastically limited here 
by air acoustic damping and a fitted adjusted effective 
mass of Meff 2 x 10“^^ kg. The nanowire can be effi¬ 
ciently driven into motion through an electrostatic actu¬ 
ation by approaching a sharp electrode polarized with a 
time-dependent voltage SV in the vicinity of the nanowire 
extremity, perpendicularly to the optical axis (see SI). 


Large deflections Sx can be achieved (see Fig. 2b), with 
a typical efficiency of 6F/6V = 50pN/V for Vq = 100 V, 
expressed as an equivalent local force applied perpendic¬ 
ularly on the nanowire extremity. For comparison, the 
nanoresonator force sensitivity (5Fth[f^] = 
amounts to 1.4fN/\/Hz in air. 

Optically resolved enhanced Brownian motion- A natural 
metric for resolving the trajectory of the oscillating sin¬ 
gle photon source is the absolute spatial resolution of the 
optical apparatus, defined via its point spread function 
n(r) (Fig. Ig). The NV fluorescence properties become 



FIG. 2. Optically resolved enhanced Brownian mo¬ 
tion. (a) Electric scheme employed for electrostatic actua¬ 
tion, combining an offset Vq and the drive signal SV. (b) 
Static deflection obtained for a ~ 10 pm electrode-nanowire 
spacing, presenting a quadratic dependence in the voltage 
applied. The force SF becomes linear in the (weak) con¬ 
trol signal when added to a large offset, (c): Employing a 
monochromatic signal resonant with the first mechanical res¬ 
onance drives a coherent oscillation of the nanowire, with a 
characteristic camelback shape visible in CCD fluorescence 
images . (d): Using a white Gaussian voltage noise thus per¬ 
mits generating an additional Langevin force which increases 
the effective vibrational temperature Tgff of the fundamen¬ 
tal mode. Scanning fluorescence images in the vertical xy 
plane for increasing excess noise strength which permits a di¬ 
rect visualization of the position distribution of the enhanced 
nanowire thermal noise (rms amplitude Aa^th)- Their horizon¬ 
tal cross sections shown in (e) are adjusted with Gaussian dis¬ 
tributions of width Wq/4: (s^^ text) and reported in 
(f) (solid lines). The green dots represent the results of direct 
noise thermometry which are only valid for spatial spreadings 
smaller than the optical waist (gray region). 
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strongly position sensitive when the spatial spreading of 
the single photon source trajectory is comparable to the 
optical waist. This regime has not yet been explored in 
existing optomechanical experiments but will soon be¬ 
come relevant with the ongoing trends towards reducing 
the oscillator mass in particular through the use of car¬ 
bon based nanoresonators and improvements of the abso¬ 
lute optical resolution through advanced super-resolution 
methods [32] . The large efficiency of the electrostatic ac¬ 
tuation, combined with the extreme force sensitivity of 
the nanowire, permits to dramatically enhance its ther¬ 
mal noise by applying an additional Langevin force SFl 
emulated by a noise generator delivering a spectrally 
white Gaussian noise (see SI). Careful piezo-positioning 
of the electrostatic tip permits to dominantly drive the 
first eigenmode in a direction perpendicular to the op¬ 
tical axis (x axis), the ID trajectory being described 
by ^(t). Increasing the r.m.s. amplitude of the noise 
source thus permits enhancing the effective temperature 
Teff of the first eigenmod e and thus the spatial spreading 
Axth = of the nanowire position fluctu¬ 

ations. This was verified at small temperature increase 
(Fig. 2f) in the direct optical readout based on real-time 
position sensing (Fig. Ig). At larger effective tempera¬ 
ture, the spatial distribution of the nanowire thermal 
noise can exceed the transverse PSF size and can thus 
directly be measured through a scanning fluorescence 
imaging (see Fig. 2d). The data are fitted by a convo¬ 
lution between the stationary Gaussian probability dis- 

tribution P{x) = (27rAx|j^) exp(—x^/2Ax|j^) associ¬ 
ated with Brownian motion and the fluorescence PSF 
of the confocal microscope sampled across the vibration 
axis. The latter which is measured in absence of drive 
(see Fig. le) is approximated at the waist by a Gaussian 
function n(x) = exp(—2x^/reo). The resulting convolved 
Gaussian spatial width /4 + Ax^j^ (Fig. 2f) allows 
measuring the spatial spreading of the nanowire thermal 
noise. The linearity of the actuation is verified up to ex¬ 
tremely large rms amplitudes, approaching the /im range. 
This corresponds to an equivalent effective temperatures 
around 10^ K for the first eigenmode and permits explor¬ 
ing nano-optomechanical systems far beyond the (ther¬ 
mal) Lamb-Dicke regime = Axth/'?^o < !)• 
Autocorrelation function- The impact of the enhanced 
nanomotion on the detected photon statistics was then 
investigated through fluorescence second order correla¬ 
tion measurement from different points in space and time, 
g^^^ (r,xi,X 2 ). This quantity reflects the probability of 
detecting a photon in the stop channel, collecting the 
fluorescence around a position X 2 , at a time difference 
r after having detected a photon in the start channel, 
centered at xi. The detection of a start event at a 
time Ti initializes the NV in its ground state and lo¬ 
calizes it at the start position near xi. Its sub¬ 

sequent fluorescence rate is proportional to the evolu¬ 
tion of the excited state population cre(r, ^(ri + r)), the 
photons being emitted from the position ^(ri + r). It 
can be numerically computed once the trajectory f{t) 
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FIG. 3. Cross-correlation function of the oscillating 
single photon source (a) measured when both detectors 
are placed at the center of the Brownian trajectory (b) for 
an increasing incoherent drive amplitude and normalized us¬ 
ing the measured photon fluxes product [33]. The strength 
of the photon bunching appearing at large excursions is re¬ 
produced in (d) and fitted with theoretical expression (see 
SI), as well as the reduction of the emitter effective lifetime, 
(c): data are well fitted with expression the fitting coeffi¬ 
cients A 2 j( 2 p 0 l )^^are shown in (e) and compared to the 
theoretical expansion (lines) converging at moderate spatial 
spreading (see SI). 


and the pump intensity profile I{x) are known. In our 
numerical simulations, a simplified 3-level scheme was 
adopted to model the NV fluorescence properties, see SI, 
with a position dependent pumping rate proportional to 
the time-varying intensity seen by the nanowire I{f{t)). 
The PSF of both measurement channels were measured 
experimentally and can be modeled by displaced Gaus¬ 
sian profiles Ili{x) = Il{x — Xi) so that the detection 
rates of the start,stop photons are weighted by ni(^(ri)), 
Il2(C(^i +^))- Here we restrict ourselves to the situation 
where the emitter lifetime (F”^) remains short compared 
to the duration of illumination, Axth/'^^o ^ r/f^rn, so 
that the probability to detect a start photon can be ap¬ 
proximated to the local time averaged fluorescence rate 
at the start position, proportional to de(^(ri)), which is 
directly measured in Fig. 2d. Finally the measured nor¬ 
malized spatio-temporal cross-correlation function of the 
vibrating single photon source, can be expressed as: 

{t,Xi,X2) = {t,Xi,X2) in, (1) 

where, using (.. .)t as the average over the in¬ 
tegration time T, we have G^^^ (r,xi,X 2 ) = 

(5-e(C(Tl))0'e(T,C(ri + T))ni (^(Ti))n2 (C(ri + 
which is integrated over all the possible start events 
(ri) for a given trajectory f{t). The denominator 
reads N = (5-e(^(t))ni(^(i)))T((fe(^(t))n2(^(i)))T and 
using the ergodic principle, we have M = 4>i4>2 with 
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= f dxF(x)ni(x) which normalizes to the product of 
the time averaged photon flux seen by each detector [33] . 
These expressions have a general reach, they are valid 
for any trajectory in space and can account for optical 
saturation or for sharp optical illumination (wq < Axth)- 
It is interesting to note that in the particular case of 
broad illumination, the excited state population recovery 
becomes position insensitive, cre(T, ^(ri + r)) ^ cre(r) 
so that the autocorrelation function can be factor¬ 
ized as: G^^\t,Xi,X 2 ) = aeir) Cislir, xi, X 2 ), with 

Gill(T,xi,X2) = f dTiIIi(^(Ti))n 2 (^(Ti + t)) Capturing 
all of the spatial dependency. This permits to clearly 
isolate the measurement-related contributions to the 
correlation function. In the following, we will explore 
different experimental configurations that permit a full 
characterization of the emitter trajectory in space. 

First, Fig. 3a represents g*^^^ (r, 0,0) measured when 
both detectors are collecting the fluorescence from the 
center of the single photon source Brownian trajectory. 
When increasing its spatial spreading beyond the optical 
resolution, a pronounced oscillatory bunching signature 
appears, whose amplitude increases with Teff. Its mag¬ 
nitude is reported in Fig. 3d and can be well fitted with 

gE (0,0,0) = / dxF(x)n(xf/^^ = (1 + 402)/yrq:^ 

(see SI). This reveals the emergence of a new photon loss 
channel for the detection path, due to photon emission 
outside of the detection volume. In analogy with photon 
bunching signatures appearing for quantum emitters 
presenting a dark metastable state, this bunching 
signature is accompanied by a sharpening of the anti¬ 
bunching profile so that the slope at origin is steepened 
from 14ns“^ to 4ns“^ at large effective temperatures. 
These signatures represent a quantitative analysis tool 
for super-resolution experiments on spatially moving 
targets. 

The demonstration of this randomization of the pho¬ 
ton emission in space due to the nanomotion was 
further investigated by measuring the fluorescence 
cross-correlations from different locations in space, 
using the movable photon counting areas. Figures 4a,b 
represent maps of (r, 0, J) (oc,0,0), where the 
start counter remains centered, the stop counter being 
displaced across the NV trajectory. Experimentally, 
it is obtained by accumulating the correlations during 
a fixed duration T at different position of the stop 
detector. The anti-bunching character of the single 
photon source is preserved all across the space direction 
while the bunching magnitude now varies in space and 
time, revealing the diffusion of the single photon source 
following its localization in = 0 at r = 0 by the start 
photon detection event. Due to the reduced quality 
factor limited by air damping, only few oscillations are 
visible before erasure of the initial conditions, when 
the autocorrelation function loses its time dependency 
and converges towards a Gaussian thermal distribution 
(see SI). These measurements permitting to visualize 
the diffusion in space and time of the single emitter are 
performed in the photon counting regime and reproduce 


the results derived from a record of ^(t) in real time on 
the QPD obtained at lower effective temperature (see 
SI). Our results are also in good agreement, see Fig. 4a,b 
with numerical simulations based on equation 0 , where 
the Brownian motion trajectory was numerically gen¬ 
erated using the experimentally determined parameters 
(effective temperature, frequency, quality factor) (see 
SI). 

Nanomotion sensing in the photon eounting regime- 
Such cross-correlation measurements in space and time 
can be employed to measure the vibration noise spectrum 
of the nanowire defined as the Fourier Transform of its 
position autocorrelation function C^{r) = {^{t)^{t + r)). 
Here we illustrate this connection between both corre¬ 
lation functions allowing to measure the thermal noise 
of the nanomechanical oscillators at ultralow optical in¬ 
tensities, falling far below the dark noise of standard 
photodiodes. Exploiting the Gaussian character of the 
nanowire thermal noise position distribution which has 
been verified experimentally (Eig. 2e), it is possible (see 



FIG. 4. Probing the oscillator dynamics with a single 
photon source. Diffusion in space and time of the oscillating 
single photon source: (a, b): the cross-correlation functions 
(t, 0,(5) (oc,0,0) are measured for different posi¬ 

tions of the stop detection area {S = 0, 0.2, 0.35, 0.5, 0.6 /xm), 
scanned across the spatial distribution of an enhanced Brow¬ 
nian motion of 0.95 ico rms spreading. Using this normaliza¬ 
tion permits revealing signal spreading in space. The results 
of numerical simulations are shown for comparison (a: full 
lines, b: bottom). Measurement of the oscillator autocorrela¬ 
tion function: c) results of the numerical simulation (see SI) 
showing the cross-correlation function g^^^ (r, S, —5) obtained 
for increasing separations ±S between the start/stop detection 
areas (see inset) and for Axth/ct;o = 0.4. d): corresponding 
experimental data obtained for similar excitation strength, re¬ 
vealing the oscillator autocorrelation function. The data are 
fitted with eq. ([^. (e): Fourier transform of g‘^{r) — 1 fitted 
with a thermal noise spectrum (dashed line). 
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SI) in the case of spatially broadened illumination to ex¬ 
pand the autocorrelation function in powers of 


{t,xi,X2) 


^0 h 


3 


( 2 ) 


The Aj coefficients depend on the spatial gradients of 
the displaced PSF and on the thermal noise spreading 
Axth (see SI). Aq reflects the product of the mean 
photon fluxes seen by each detector and thus con¬ 
veys the role of the shot noise in that measurement. 
Maximizing Ai with respect to higher orders terms 
allows to optimize the measurement of the nanowire 
autocorrelation function C^(t), which is not possible 
when one detector remains centered (such as in Fig. 
3,4a,b) since ^ 2 j+i = 0 there. Instead, it is interesting 
to reproduce a quadrant photodiode configuration, 
when both photon counters are arranged to monitor the 
fluorescence on each side (±6) of the nanowire trajectory 
(see Fig. 4c). In that situation, larger contrasts can be 
obtained for sufficient separation (see Fig. 4c), and the 
quadratic term A 2 can even be nulled for 6'^ = wo/2 
(see SI). The NV cross-correlation function now presents 
an oscillating and exponentially decaying pattern, 

see Fig. 4d, which can be perfectly adjusted using 
equation ^ and the position autocorrelation function 

C^{r) = , 

with — r^/2 [34]. Its Fourier transform 

(Fig. 4e) provides a measurement of the enhanced 
thermal noise spectrum Sx[^] of the nanowire. The 
emitter finite lifetime is responsible for an upper bound 
of ^ 40 MHz on the measurement bandwidth which 
does not limit our measurements here. For small 

excitation amplitudes (Axth one obtains: 

(5*, —(5*) o-e{r) (1 — /wq) ^ (see SI) with a 

shot noise limited sensitivity of Wq/ (44>yT^) amount¬ 
ing to (6nm)^/\/Hz for a 10^ Hz mean photon flux and 
a bin time of Tbin = 1 /^s. This sensitivity is sufficient 
in principle to detect the thermal noise of ultralight 
oscillators such as suspended carbon nanotubes, even at 
low temperatures (see SI) [35] . 

We note here that it is generally possible to realize a 
linear measurement of the oscillator position ^(t) by 
recording and subtracting the photon fluxes seen by 
each detector arranged in a quadrant photodiode config¬ 
uration. In that case, the thermal noise can be detected 
above the photon shot noise if 4> > Tjji{wo/Axth)‘^ which 
requires having large photon fluxes, light oscillators or 


high mechanical quality factors. However this method 
suffers important experimental constraints since in order 
to convert the photon counting signals into a meaningful 
time-resolved intensity signal, it is necessary to choose 
an integration time comparable to the inverse mean 
photon rate (4>“^). In order to probe the nanomotion in 
real time, this integration time should be smaller than 
the oscillation period, so that it requires in turn to detect 
at least one photon per oscillation period In 

contrast, our experiment operates below this criterion, 
which underlines the strength of this method based on 
second order fluorescence cross-correlations, even for 
detecting high frequency oscillators. 

Conclusions- This proof-of-principle experiment 
demonstrates the possibility to detect and analyze 
the vibrations of nanomechanical oscillators in the 
photon counting regime by recording the spatial and 
temporal correlation functions of the emitted photon 
flux. The detected optical powers involved fall in the 
sub-fW range, where standard analog detectors are 
dark-noise-limited. Furthermore this permits to operate 
at photon fluxes smaller than the oscillation frequency, 
where time-resolved position measurements are not 
possible with analog sensors. Therefore our approach 
also proves useful when very limited optical powers 
must be employed, such as in cryogenic experiments. 
The exposed measurements can virtually be transposed 
to any nanomechanical system that is weakly coupled 
to light (such as nanotubes) by working on Raman 
scattered photons [36] or on other defect-related fluo¬ 
rescence signals |36ll38] but also to optically trapped 
emitters [SD ESJ. The sensitivity of the measurement 
can be significantly improved through the use of optical 
super-resolution techniques [32| or microwave assisted 
sharpening of the spin fluorescence PSF gniST]. The 
investigation of the qubit dynamics second order corre¬ 
lation function can also be extended to the frequency 
domain of parametrically coupled hybrid mechanical 
systems, through spectral diffusion analysis [42], or 
resonant optical pumping [a EH na sa, where larger 
coupling strengths can be achieved. 
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Supplementary material for “Nano-optomechanical measurement in the 

photon counting regime” 


Experimental setup 


White 



Spectrum Analyzer 


Computer Acquisition 


Gating and Counting 


FIG. SI. Simplified scheme of the experimental setup. BS: Beam Splitter, QPD: Quadrant Photodiode, APD: Avalanche 
PhotoDiode (single photon counting modules). 


This paragraph gives a more detailed description of the setup shown in|Sl[ The laser is a Laser Quantum “Gem” 
laser (532 nm). The optical power, size and polarization of the beam can be adjusted prior to its injection in the 
microscope objective. The optical power is controlled with a combination of a multi-order half-wave plate and a 
polarizing beam-splitter. A telescope is used to match the beam diameter with the size of the input iris. Another 
half-wave plate is used for polarization control. 

The beam is then focused in a long working-distance (4 mm) and high numerical aperture Zeiss microscope objective 
(xlOO, NA=0.75). After interaction with the functionalized nanowire, some of the back-scattered light is collected 
through the same objective. 

A Semrock 580 nm dichroic mirror serves to separate the intense green excitation beam (about 0.5 to 3mW) from 
the weak red fluorescence outcome (about 1 fW) from the sample. This dichroic mirror stands just in front of the 
microscope objective and reflects the green laser into the objective while it lets the red fluorescence pass through and 
reach the fluorescence analysis stage. The backscattered green light is also sent back via the same optical path as 
the injected laser to be analysed. A third collection path is set up by an identical microscope objective positioned 
opposite the injection objective to collect the transmitted green light or the fluorescence. 

A white light source is also injected through the backside of one of the injection mirrors. Its reflection on the sample is 
collected through the dichroic mirror and directed with a mirror to a Watec 910HX camera which permits fluorescence 
imaging and sample monitoring. 

Additional filters are added after the dichroic mirror in order to completely suppress the pump light in the fluorescence 
analysis path. It is split with a non-polarizing beam-splitter into two ideally equally intense optical paths fed into 
two Avalanche PhotoDiodes (Single Photon Counting Modules) through two 50 jam core multimode fibers. The 
fluorescence signal is converted into voltage pulses by the APDs which are in turn either gated and counted by the 
interface program to perform nanowire positioning and tracking or fed to a FastComTec MCS6A correlator to build 
intensity correlations. One of the fiber couplers is mounted onto a PI 500 jam range lateral piezo axis so as to move it 
in the image plane of the sample and to build intensity correlations from different points in space. This supplementary 
axis was first used to determine the real magnification factor of the experiment ( [S^ . 

An asymmetric (90:10) beam splitter reflects most of the reflected light towards a Quadrant-PhotoDiode (QPD) 
with a homemade differential amplifier. The difference and the sum of the voltages of the two quadrants of the 
photodiode are measured and both the DC and FIF components are split, thus resulting in four output channels. The 
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FIG. S2. Calibration of the real magnification factor. Fluorescence map obtained by scanning the object and movable detector 
positions. The laser light field is broadened by a defocussing lens placed before the objective, which enables the NV center 
to be lit over several microns. The solid black htted line yields a magnihcation factor of 128.2, in good agreement with the 
expected value 125 for the combination of focalization lens and objective used here. 


HF component is fed into a Spectrum Analyzer (Agilent MXA) to build the position correlation of the nanowire and 
the DC signal is fed to the interface program for imaging purposes. 

The mechanical oscillator is a 46 nm-long silicon carbide nanowire attached to an electrochemically etched tungsten 
tip. The etching, preliminary characterization and binding of the nanowire to the tip are performed in the ILM 
in Lyon in collaboration with A. Siria, P. Vincent, A. Ayari and P. Poncharal. A XYZ Piezostage supporting the 
sample holder of the nanowire can be moved in front of the microscope objectives so as to scan maps of the reflected, 
transmitted and fluorescence lights as a function of the nanowire position in the optical beam. Generally, two axes 
are scanned to produce ”XY” 30 x 30 jam or ”XZ” 30 x 10 jam maps. These scans and the analysis of all signals are 
performed thanks to a NI e-6323 measurement card and a homemade Python-Qt interface adapted to the setup. To 
get the laser waist to lie in the 30 x 30 x 10 fim piezo window, the whole piezostage is mounted on a XYZ Newport 
translation. 

Nanodiamonds are attached to the nanowire by dipping the nanowire into a droplet of commercial solution containing 
50 to 100 nm diameter nanodiamonds. They are small enough to contain 0 or 1 defect at most, given the typical NV 
defect densities of this type of diamond. 

Mechanical vibrations 


The vibrations of the nanowire are probed with a focused laser beam illuminating the NV defect, which is located 
at the nanowire extremity. The description of the optical readout of the deformations of the nanowire is presented 
in m, where the multi-modal expansion is presented in detail. It permits defining an effective mass, see |S2j . which 
depends on the vibrational mode considered and on the position along the nanowire. In case of a small optical waist 
with respect to the length L of the nanowire (and to the other characteristic sizes of the mode deformation profile) 
the optical readout permits a local measurement of the nanowire deflection. The effective mass can be expressed as: 


Meff,n = at 



dy Un^{y) 

L uJiL)^ 


(SI) 


where Un{y) is the displacement of the beam at longitudinal position y. y = L corresponds to the free extremity and 
y = 0 the clamped extremity of the nanowire. The numerical resolution of the Euler-Bernoulli equation for a singly 
clamped beam gives the non-normalized shape of the spatial modes of the nanowire (see Table [ij), which in turn allows 
to calculate the effective mass of each vibrational mode (see Table |^. We have: 


(t) = cos{kny) - cosh{kny) 




sm{kny) - smh{kny) 


}■ 


(S2) 
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Position (ijm) Excitation ampi. (mVpp) 


FIG. S3. CCD fluorescence images obtained under broadened spatial illumination for increasing coherent driving of the 
nanowire. The Airy peak (left) observed at rest is progressively broadened, presenting a characteristic camel back shape. 
The driven trajectory can spread over several fim. To do so a resonant RF tone is applied on the electrostatic tip, with 
amplitude from 0 to 45 V (after the voltage amplifier) and offset of 20V. The prohles are htted with an oscillating Gaussian 
function: F(x) — ^ ^- 2 (x-xo-Axcosnmt) averaged over one mechanical oscillation (period Tm) centered in xo with 

an amplitude Ax. These two quantities are reported for increasing drive amplitudes and fitted with quadratic an linear functions 
respectively. 


n 

kfiL 

A-n, 

Meff/M 

1 

1.87510 

-0.7341 

0.2500 

2 

4.69409 

-1.0185 

0.2500 

3 

7.85476 

-0.9992 

0.2433 

4 

10.9955 

-1.0000 

0.9547 

5 

14.1372 

-1.0000 

0.9646 


TABLE I. Numerical coefficients corresponding to the first 5 eigenmodes of equation |S2| in case of a point-like optical mea¬ 
surement (homogeneous mode deformation profile within the waist area). 


This work has been realized at ambient pressure, where the two orthogonal polarizations of each longitudinal eigen¬ 
mode family are not resolved due to air damping. In addition, the nanowire is always positioned on the optical axis, 
where the quadrant photodiodes are only sensitive to the nano-motion perpendicular to the optical axis. Furthermore 
the electrode is micro-positioned perpendicularly to the optical axis, so that it generates a force that is driving the 
nanowire perpendicularly to the optical axis. All these reasons justify the simplification of the description of the 
nanowire of the nanowire vibrations by a uniaxial oscillator, only oscillating perpendicularly to the optical axis. A full 
description of the nanowire dynamics in 3D can however be found in [ST], but is not required here to fully describe 
our results. In the following and in the article, we then describe the measured transverse position of the nanowire by 
the scalar Sx{t). 

Furthermore, we will only restrict our study to the fundamental eigenmode. This is fully justified since the second 
longitudinal eigenmode oscillates around 6 times faster than the fundamental period. This means that compared to 
the fundamental mode, the r.m.s. Brownian motion spreading of the second mode at the extremity is around 36 
times smaller and the noise power of thermal noise (oc Q/M/ft^) 216 times smaller at resonance. Moreover, the 



















4 


spatial profile of the electrostatic actuation is dominantly exerted at the extremity of the nanowire, in front of the 
electrostatic tip, so that the relevant increase of the effective temperature of the higher order modes is significantly 
reduced as compared to its effect on the fundamental eigenmode. This is even more pronounced due to the finite 
bandwidth of the noise generator that can be tuned to fade at higher frequencies, without impacting on the white 
character of the noise seen by the fundamental mode. The mechanical susceptibility x[^] is defined as: 

Sx[Q] (S3) 


where is the spectral component of the displacement measured at the free extremity of the nanowire and 

is the spectral component of an external force. Following the previously mentioned approximations, the nanowire 
dynamics can be assimilated to a single mode oscillator with the susceptibility: 


yfril — 

(S4) 

The mechanical displacement spectrum is defined as: 


s,[n] = \xm^sF[n], 

(S5) 

with 


Sf[^] = 2MeffrmA^BT, 

(S6) 


where T is the effective vibrational temperature of the nanowire. 


The position of the nanowire Sx{t) can be read out in real time and recorded through the reflection or the trans¬ 
mission of the green laser light onto one of the two QPDs of the setup. It is then possible to compute the oscillator 
space-time correlations Csx{^) = {Sx{t)6x{t + r)) as shown in Fig. 


S4 


Since the readout laser is positioned on the NV location, the measured displacement Sx{t) also represents the transverse 
position of the NV defect ^(t). Once again, we restrict our analysis to a mono-dimensional trajectory, perpendicular 
to the optical axis. The extrapolation to any other 3D spatial trajectory is straightforward once the PSFs are 

known. 
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FIG. S4. Space-time correlations of the position of the nanowire extremity undergoing Brownian motion at 300 K, measured 
using the QPD detectors (green trace). To compute this correlation a long lasting signal (about 3 s) from the QPD measurement 
apparatus is acquired (12-bits resolution) and calibrated in displacement. 1800 sections of 30 fis are singled out if they meet 
the two following criteria: start at a position around 3.8 nm from the equilibrium position and with a speed around 0nm.s“^, 
with both selection windows set at the limit of resolution of the apparatus roughly given by the trace initial breadth. These 
criteria set a start point in phase space from where we observe its re-population under decoherence. The chosen criteria are 
therefore arbitrary and the figure serves only as an illustration for this measurement [S3] . Note that an additional condition 
is added for the computation which is that no subsequent sections separated by less than 30 fis are retained in order to avoid 
parasitic correlations due to over-sampling. Eventually each of the chosen sections is digitized over a grid, and the sum of these 
grids give the correlation trace presented here in green, thus rebuilding in time domain part of the signal usually computed by 
a spectrum analyzer. The space average position (dark red solid line) is computed from the analysis of vertical slices of this 
trace. 

Alternatively, the high frequency output of the QPD is sent to a Spectrum Analyzer (Agilent MXA) to compute 
the position spectrum. The final position spectrum showing several vibration eigenmodes of the nanowire in Fig. Ig 
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was extracted from the Spectrum Analyzer raw measurement 5 'vhf [^] through: 

where 5'darkHF the (dark) background noise of the QPD HF channel, Gdc[^] and Ghf[^] are the previously 
calibrated DC and HF frequency-dependent gains, and is the DC-voltage-position conversion factor measured 

by making a DC image of the nanowire displaced with the piezoelectric stage. Typically, the laser shot noise emerges 
above the QPD dark noise level for optical powers larger than ^ 100 — 200 /iW. 

Electrostatic actuation 

A conducting tungsten tip is piezo-positioned in the vicinity of the nanowire’s extremity. A high voltage - 50- 
fold amplification of an Agilent Arbitrary Waveform Generator typically 3V signal - is applied between the tip 
and the nanowire metallic holder. Other than this tip the whole setup is electrically connected to the holder and 
grounded. The electrostatic force is the spatial gradient of the energy stored in the circuit formed by the nanowire 
and the electrostatic tip. For a tip positionned perpendicularly to the optical axis, the expression of the equivalent 
electrostatic force applied on the apex of the nanowire can be written as: 

F(ro) = -aV^e^, (S8) 

where V is the applied voltage difference and a contains all of the spatial and material dependence of the electrostatic 
description of the system. 

According to classical beam theory, the static deflection of the nanowire is then in turn quadratic in voltage as 
shown on Fig. 2b: 


\5x\ 


aV'^F 

3EI 


kV^ 


(S9) 


where L is the nanowire’s length, E its Young modulus, / its moment of inertia, and is a global electromechanical 
coefficient for the system. 

The fit of the static deflection as a quadratic function of voltage shown in Fig. 2b gives an electromechanical 
coefficient of /^stat = 8-53 • 10“^^ m.V“^. 

The setup illustrated in Fig. 2a is used to simulate a Langevin force of arbitrary temperature with a white Gaussian 
voltage noise SV of spectral density and variance ay added to a voltage offset Vq. This offset is chosen to be 

much greater than the force signal so that at first-order the dynamic force noise becomes linear in the noise voltage 
(at the cost of a high static force and thus a strong static deflection which is compensated by a displacement of the 
nanowire support): 6F ^ 2aV{)5V. The corresponding force spectral density becomes 




and the effective temperature (the damping rate is unchanged): 


Teff = T + 


2a^V^ Sy 


(SIO) 


(Sll) 


The signal generator bandwidth is chosen to fall short of the second mode frequency. Only the first eigenmode of the 
oscillator is sensible to this simulated high temperature. The flatness of the generated voltage signal was verified with 
a Spectrum Analyzer and its Gaussian character was also assessed through quadrature analysis for several frequency 
ranges. 

The time-averaged spatial distribution of the NV center fluorescence <F(xo) is measured with an APD (Fig. 2a) 
under wide illumination field. When the nanowire is driven with an enhanced and spectrally white force, the r.m.s. 
position can be inferred from a Gaussian fit of the time-integrated fluorescence : 


^{xo) = 


/ + 00 
-oo 


n(x — xo)P{x)dx = <F 


/ + 00 

' 

-CX) 


2 (x-xq) 

> 


th 


<Foe 

^l + AUxlJwl 


--dx 


(S12) 


where wq is the detection waist and Axth the r.m.s. position. 
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FIG. S5. Initial bunching/anti-bunching strength Golc(0, xi, a: 2 )/Golc(oo, xi, X 2 ) of motional origin obtained for varying 
start/stop detector positions ^ 1,2 = xi, 2 /'?r’o/v^ in the case of increasing oscillating amplitude Ox = 0.1, 0.5,1.5 (left to right). 


The orange dots in Fig. 2f show the total fitted Gaussian width of the broadened fluorescence images. The red 
line shows the same width corrected for the detection waist to give Axth- This is therefore a second method for 
calculating the electromechanical coefficient : K^yn = 6.70 • 10“^^ m.V“^ which coincides reasonably well with the 
static characterisation. Yet another measurement of the r.m.s. position is given by the integration of the position 
spectral density (green dots on Figure 2f), which is in agreement with the previous measurement as shown by the 
equal slopes between red and green dots at small effective temperatures on Fig. 2f. However in this measurement 
the position was measured from the QPD differential voltage which has a limited linearity zone that the nanowire’s 
motion exceeds when the temperature is increased, resulting in saturation effects. 

Generic expressions of the initial autocorrelation strength 

Using the spatial distribution of the Brownian motion 

1 

P(x) = , -.e 2Aa;2 

V‘27rAx‘^ 

and the displaced PSF,: 

ni(x) 


we can estimate the initial bunching. We have: 

G^S{oo,Xi,X2) oc /<feP(x)n.(x)/<!xP(x)n,(x) 


and 


G'is^( 0 ,a;i,a; 2 ) oc / dxP{x)Ili {x)Il 2 (x) 
so the initial relative correlation strength gives, using Si = Xilw{^l\f2\ 


Gosc(0,xi,X 2 ) _ 1 + 4^^ 


G'is^(oo,a;i,a;2) \/l + 80® 


exp 


ml 


(1 + 402)(i + 802) 


(202(5^ _ ^2)2 _ 5^52) 


It is plotted in S5, for varying vibration strength. In the centered case one obtains: 

dskoy^o) _ i+4^g ^ ^ 

gE(oo, 0,0) yi + 802 

which is in good agreement with experimental results (see Fig. 3d of the manuscript). 


(S13) 


Simulations 
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As spin variables here are irrelevant, the NV center is modeled as a three-level system, keeping only one ground 
state 1^) and one excited state |e) regardless of spin projection, and merging the whole S = 0 levels system 
into one dark metastable state \m). The statistical population of these states will be denoted de and 
and normalized such that the sum of these three populations is always unity. The fluorescence intensity is then 
proportional to the excited state population ag. Right after the emission of one photon, the state of the NV center 
is known to be \g). Thus, the intensity correlation can be computed from the evolution of (jg under the optical 
Bloch equations with the initial condition of ag = l,crg^^ = 0 as crg(t) then represents the probability of emitting 
a second photon. Any (slowly) time-dependent excitation intensity can be injected in the Bloch equations for 
this numerical resolution which allows us to build the intensity correlations for any type of trajectory of the NV 
center in and out of the excitation volume. Additional spatial weighting is added to model the motion in and out 
of the start/stop collections volumes which may or may not coincide, depending on the type of measurement performed. 

In adiabatic cases where the NV center’s populations adapt almost immediately to the received optical power 
compared to the mechanical oscillation period, the resolution of the Bloch equations is not essential and the 
fluorescence, even modulated by the slow motion of the nanowire, can be computed either analytically for simple 
trajectories as sine-shaped ones or through statistical considerations for Brownian motion. However, in order to 
take into account the photodynamics of the NV center and investigate the non-adiabatic regimes using a simulation 
Runge-Kutta 4 algorithm to resolve optical Bloch equations. 


Since the choice of trajectory is then completely arbitrary, this program can be fed with a simulated trajectory 
obtained with the convolution of a white noise are then used to build the full correlation. 

Each of these simulated trajectory realizations is first given a weight depending on the presence or not of its starting 
point within a Gaussian window around the first detector’s central position. This accounts for the start channel PSF. 
The subsequent probability of fluorescence through time is calculated by numerically solving the Bloch equation for 1 
point out of 2 in the trajectory (as the Runge Kutta algorithm requires twice as many points as wanted in the output 
signal). Each of these fluorescence values is in turn pondered with a second weight to take into account the Gaussian 
window around the second detector’s central position. The two detectors can therefore be moved separately. 

Expansion of the photon autocorrelation 


We show here that in case of broad illumination, the photon autocorrelation function can be expanded as: 


/ \ oo 


10 


j=0 


wll2 


(S14) 


where the Aj coefficients are delay independent and functions of the detector positions Xi and oscillator spatial 
spreading Ox- One has: 


A. 


2i 


(2i) 


E 

n,m=0 


{2n + 2j)\{2m + 2j)\ 
{n +j)\{m +j)\ 


-2xi 


-2x1 


m 

A even 


and 


^2i+l 


-2xiX2 ( 2n + 2i + 2 )!( 2TO + 2j + 2 )! /- 2a;fy dd 

wU2j + l)\^^^^ (n+j + mm+j + l)\ \ J \ J 


with 


A even 


E 


i^Y (^) 


(2n - 2 p)!(2to - 2q)\p\q\ ’ 


a: 


odd 


= E 

p,q=0 


[2xj J 2a| ) 


{2n -2p + l)!(2m -2q + l)\p\q\' 



In case of broad spatial illumination, the measured fluorescence is modulated on short time scales by the photo¬ 
physics of the NV center and on long time scales by the motion of the NV center in and out of the collection volumes. 
Rewriting the expression introduced in the manuscript in case of broad illumination gives: 


{t,Xi,X2) 


<^e{T) (n(^(t) - Xi)n{^{t + r) - X2))t 

N 


(S15) 


We will now expand in powers of ^(t) the two probabilities of detection centered respectively on xi and X 2 , and 
in doing so we will only suppose that the PSFs are expandable in power series. We will focus on the numerator 
G^‘^\t,Xi,X2). 


(t,Xi,X 2) = (Je(j) ( y] y] 


(T-n 


\n=0 m=0 


dx^ 


d^U 


dx^ 


{^{t) - Xi)^{C{t + t) -X 2 Y' 
nlml 
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{t,Xi,X 2) = (Teir) ^ ^ 


1 d^U 


n=0 m=0 


nlml dx^ 


d^U 


dx'^ 




*%=0 9=0 \P 



^-^xrHatnit+ry) 


(S17) 


where we have omitted the subscript T in the time averaging. 


denotes the m-combination from an ensemble 


of size n. Since ^(t) is a Gaussian distributed random variable, as visualized in Fig. 2 of the manuscript, a very 
useful property of Gaussian processes can be used to reduce -h tY) to a sum of small manageable terms. 

This property holds as follows: for n jointly Gaussian random variables a::{i=o...n} whatever their auto- and eross- 
eorrelations: 


(Xl-Xn) = 
{xi...Xn) = 0 


if n is even 
if n is odd 


(SI 8 ) 


where Cpairs denotes the ensemble of pair configurations for the indexes l...n, where each index appears exactly one 
time and with no consideration of order. 

This property therefore holds for our random variables ..., ^{t)^ ^{t-\-T )^..., ^(t+r)} and enables the development 


of {^^{t)^^{t + r)) into terms of the form: C|(r)Ax^^. In order to amount the contribution of each term to the total 


correlation, a diagrammatic picture of the problem was built upon the principles illustrated in S6 In this picture, a 


vertex represents the position at a time t, ^(t), and a line joining two vertices represents a correlation between these 
two positions, that is, autocorrelations of the oscillator position. A line can also join one vertex to itself making a 
loop which thus represents the zero-delay autocorrelation. 


• 

• 

\b) 1 

Cdr) 



i{t + T) 

^{t) ({t + r) 

0 


FIG. S6. Building blocks of the diagrammatic representation of terms, a) Dots are positions of the oscillator at different times, 
b) Lines that connect one dot to another represent autocorrelations, c) Loops connecting one dot are zero-delay autocorrelation. 


Some examples of this representation are given in S7 for the first two orders in correlation of the expansion. 
Gounting all the possible ways to obtain a term where the correlation appears at a fixed power j : 
from the term {^'^{t)^^{t + r)) with no consideration of order, we obtain that the possibilities are numbered as 
explained in[S8{ 


j 7 ! (g-i)' 

W W ' 

I II III IV 


(S19) 


I is the number of ways to select j elements from each ensemble of size p and q of variables ^(t) and + r) to 
produce a hybrid pair Once j elements from each group are selected, there are still jl ways to combine them to 
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FIG. S7. Example of terms reprsented with the diagrammatic picture explained above, a) Non exhaustive list of diagrams 
contributing to the first power of correlation in the total expansion, b) Non exhaustive list of diagrams contributing to the 
second power of correlation in the total expansion. Each of these diagram is weighed by a factor calculating in the following. 



S8. Diagrammatic explanation for the calc ulation of the degeneracy factor of each diagram representing a term 
(p+g)- 


-2j 


(complementary to equation S19). Starting from the knowledge that each vertex must be connected p 


EIG. 

C'|(r)Aa:J^ _ 

(q) times and that j lines must connect them, there is a number I of possibilities to select j branches from each vertex to build 
interconnections. The lines are colored to identify them when it is relevant. After this selection of j branches per vertex, there 
is still a number II matching possibilities. Then the remaining branches are connected to other branches from the same vertex, 
which makes up another degeneracy factor III for one vertex and IV for the other. 


make which is represented by the factor IT Identically, the remaining p — j and q — j elements in each group can 

be associated two by two in different ways which results in degeneracy factor III and IV. Note that j, p and q have the 
same parity: the terms that would not respect this condition do not exist in the development. After simplifications 
one finally obtains for p and q of the same parity: 


{atfiit + rY) 


min(p,g) 

j: 

j=0 


p\q\ 




j of p,q parity 


(S20) 


and + r)^) = 0 if p and q have different parities. 

Now assuming the collection volumes IIi and 112 to be Gaussian (more precisely, TEMOO profiles at waist), centered 
in xi and X 2 , and normalized so that the maximum detection is unity. 


ni(:r) = 


2(x-xi)'^ 


^ 2 {x) = e 


2(x-X2)'' 


(S21) 


one can finally compute from equation (S17), using a serial expansion of , the total development of the numerator 
of the autocorrelation function. Finally, it is of special interest to rewrite this deve lopment as a sum over the powers 
of Ci , which results after rearranging the terms in the numerator of equation 


S14 


Once the numerator is expanded, we can point out that at infinite delay r, the mechanical autoccorelation function 
converges towards 0, so that we are only left with the Aq term of the series expansion. As a consequence, the 
denominator equals Aq. 

Let us mention here that very efficient algorithms based on cumulants can be exploited to directly compute the 
overall oscillatory autocorrelation function G^lc (r,xi,X 2 ). [S4] . 


Symmetric case (r, 5, —5) 
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F IG. S 9 . Aj expansion in the symmetric case (t, (5, —(5) a): Strength of the Aj terms calculated from the series expansion 
((S22), (S23)), as a function of the normalized temperature Ox and detector position S showing the cancelation of A 2 in 5 = 1 /\/2 
at small excitation, b) cross section for Ox = 0.001 c) evolution of AjjAo for S = lly/2 as a function of This validates the 
approximated expansion and at low drive excitation. 


We study here the case where both detectors are symmetrically positioned, monitoring each side of the trajectory: 
xi = —X 2 = We introduce the normalized detector position and oscillator temperature: 




wqI^/2 


and Ox 


Ax 

Wo ■ 


One then obtains: 


and 


A2j{6,-6) 


1 

m 



(2n + 2j)! 
{n + j)\ 


(2n - 2p)!p! 


n 


Aij+i{5, -5) 


(2n + 2j + 2)! {^x/^ ) \ 

(2j + l)! (n + j + 1)! V ) ^ (2n-2p+l)!p!j 


(S22) 


(S23) 


For small oscillation amplitudes, 0 <C 1, one can restrict the second finite sum to the power p = 0. The Aj coefficients 
are then temperature independent and one has: 
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1 / ~ (2n + 2 - 1) (-^2^ 


^ \ n \ 2 


= - 4 


i n=0 


d{-6^) 


{-S^e-^ ) - 2e-" = 2(1 - 2SYe 


^ 2 \ 2 - 26 ^ 


From that it is clear that A 2 is suppressed for 6 = 6'^ = i-e. 5 = wq/2. One then obtains: 

~ 1/e , Ai((5", -(5") « 2/e 


For higher spatial spreading Ox^ the previous simplification is not possible. The dependency of the Aj coefficient in 5 
and Ox is shown in[^ It is still possible to cancel the A 2 term by choosing the proper detector displacement. The 
numerical quantities have been calculated up to Ox = 0.35, beyond which numerical convergence is not achieved. 


At low excitation, one then obtains: 


9^'^\r,5*,-S*) «cre{T) (1 -4^1^ + ../ . 

V '^0 / 

The mechanical spectrum is then convoluted with the response function corresponding to the Fourier transform of 
cre(r), which behaves as a low pass filter of cutoff frequency (1 /tc). The number of detected photons in each time bin 
after an integration time T is ^i^ 2 ^binT, where Tbin is the bin time. For a Poissonian distribution, assuming identical 
fluxes on each APD, the equivalent sensitivity of the apparatus is then given by 

for a photon fiux of 10^ cps (300 fW at 700 nm) and a time bin of Tbin = 1 A^s. This value can be compared favourably 
to the thermal noise of a 1 MHz (Meff = 10“^^ kg) suspended carbon nanotube: Axth = /Mfl^ « 102 nm at 

300 K and 12 nm at 4K. In the article, all the measurements have been carried out at a 1.6 ns bin time to resolve 
the anti-bunching signatures - at a price of a longer acquisition time (1 hour for Fig. 4d) - which is not necessary to 
resolve the oscillating bunching pattern of motional origin. 


Comparison between simulations and Aj expansion 


Here we verify that the expansion in powers of the oscillator autocorrelation function C^{r) gives results in agreement 
with the direct numerical simulation of the oscillating NV fluorescence properties, introduced above. The comparison 
is presented here for a given excitation Axth/'?^o = 0-3 and for increasing distance between the detection areas, in the 
symmetric case: xi = —X 2 = S. The results of numerical simul ation a re sho wn in |SlQ| where they are also compared to 
the expansion up to the 4^^ order in C^(t), using expressions (|S22|), (S23) and similar photophysical and mechanical 
parameters as the one employed in the numerical simulation. The numerical simulations are also fitted with: 


(r, S,-S) 



(S24) 


SIO 


were a very good 


the only fitting parameters being the aj. They are compared to the coefficients Aj 

agreement is found. The slight deviation visible in the temporal traces obtained at large detector separation 6 might 
be due to the rarefaction of the detection events which requires larger computational strength, and to the fact that 
the expansion to the 4th order may become limited for perfectly describing this situation. 


Centered case: 0,0) 

Notice that ^ 2^+1 ^ 0 when Xi 0. For even coefficients, since 6 appears at the power the only remaining 

coefficient in the second sum is p = n which permits to express: 


^2, (0,0) 


1 (2n + 2j)\ fn 2 \ 

(2j)! (n +j)!n! 
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FIG. SIO. Comparison of the numerical simulations to the Aj expansion. Left: numerical simulations of (t,S,—S) for 
= Axth/'w;o = 0.3 (gray) for varyi ng de tector positions <5 = 0, 0.64,1.95 from i to Hi The dat a are fitted with an expansion 
up to the 4^^ order in c(t) (equation ( S24 )) (red dashed traces). Also shown in blue is expression (S14) where the Aj coefficients 
were computed using expressions ( |S22|S23 ). Right: comparison of the fitted aj(S, —S)/ao(0, 0) coefficients (connected symbols) 
and the respective coefficients Aj(S, —S)(2^x)~^ (dashed lines), which have been numerically evaluated using the 500 first terms 
of the series. The good agreement achieved validates the expansion (S14), even for strong driving. 


For simplicity, data were fitted with a normalized autocorrelation function: c(r) = C^(t)/Ax^: 

sixp (t". i -^) = ^e(T) — aj c(Ty . 

The connection between both coefficients is then: 


j=0 




2\J 


Space-Time diffusion case: (r, 0, 5) 

In that case, using the same normalized quantities and following the same reasoning gives: 


A in ^ (W + 2j)! (2n + 2j)\ 

A 2 j{ 0 ,S) (2j)! (n + j)!n! ^ J (n + j)\ \ 


n^O + 


^ {2n-2p)\p\ 


and A 2 j+i( 0 , ^) = 0 

In both situations (r, 0, 0) and g^^^ (r, 0, 5)) the nullity of odd coefficients ^ 2 j+i renders the motion dependent 

(2) 

part of the photon autocorrelation functions ^osc(^, 0,0) an even function of C^{r). Thus it does not oscillate above 
and below its limit at infinite delay. This is particularly visible in Fig. 3 and 4a and also gives rise to a contrast 
inversion which is visible in Fig. 4a, 4c. 
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